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Abstract 

We investigate inflation driven by the evolution of highly excited quantum states within 
the framework of out of equilibrium field dynamics. These states are characterized by a non- 
perturbatively large number of quanta in a band of momenta but with vanishing expectation 
value of the scalar field. They represent the situation in which initially a non-perturbatively large 
energy density is localized in a band of high energy quantum modes and are coined tsunami- 
waves. The self-consistent evolution of this quantum state and the scale factor is studied ana- 
lytically and numerically. It is shown that the time evolution of these quantum states lead to 
two consecutive stages of inflation under conditions that are the quantum analogue of slow-roll. 
The evolution of the scale factor during the first stage has new features that are characteristic 
of the quantum state. During this initial stage the quantum fluctuations in the highly excited 
band build up an effective homogeneous condensate with a non-perturbatively large amplitude 
as a consequence of the large number of quanta. The second stage of inflation is similar to the 
usual classical chaotic scenario but driven by this effective condensate. The excited quantum 
modes are already superhorizon in the first stage and do not affect the power spectrum of scalar 
perturbations. Thus, this tsunami quantum state provides a field theoretical justification for 
chaotic scenarios driven by a classical homogeneous scalar field of large amplitude. 



1 Introduction 



A wealth of observational evidence from the temperature anisotropies in the cosmic microwave 
background favor inflation as the mechanism to produce the primordial density perturbations]^, |2[. 
Thus, inflationary cosmology emerges as the leading theoretical framework to explain not only 
the long-standing shortcomings of standard big bang cosmology but also to provide a testable 
paradigm for structure formation The recent explosion in the quantity and quality of data 

on temperature anisotropies elevates inflation to the realm of an experimentally testable scenario 
that leads to robust predictions that withstand detailed scrutiny ||l], |2[. 

However at the level of implementation of an inflationary proposal, the situation is much less 
satisfactory. There are very many different models for inflation motivated by particle physics and 
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most if not all of them invoke one or several scalar fields, the inflaton(s), whose dynamical evolution 
in a scalar potential leads to an inflationary epochQ-Q. The inflaton field is a scalar field that 
provides an effective description for the fields in the grand unified theories. Furthermore there is 
the tantalizing prospect of learning some aspects of the inflationary potential (at least the part 
of the potential associated with the last few e-folds) through the temperature anisotropics of the 
cosmic microwave background 

Most treatments of inflation study the evolution of the inflaton via the classical equations of 
motion in the scalar potential and the effect of quantum fluctuations is typically neglected in the 
dynamics of the inflaton. Furthermore, since inflation redshifts inhomogeneities very fast, the 
classical evolution is studied in terms of a homogeneous classical scalar field. The quantum field 
theory interpretation is that this classical, homogeneous field configuration is the expectation value 
of a quantum field operator in a translational invariant quantum state. While the evolution of 
this coherent field configuration (the expectation value or order parameter) is studied via classical 
equations of motion, quantum fiuctuations of the scalar field around this expectation value are 
treated perturbatively and are interpreted as the seeds for scalar density perturbations of the 
metric ||-@. 

A fairly broad catalog of inflationary models based on scalar field dynamics labels these either as 
'small field' or 'large field' 0. In the 'small field' category the scalar field begins its evolution with 
an initial value very near the origin of the scalar potential and rolls down towards larger values, an 
example is new inflation||3|, ^. In the 'large field' category, the scalar field begins very high up in 
the potential hill and rolls down towards smaller values, an example is chaotic infiation|^, ^. 

It is only recently that the infiuence of quantum fluctuations of the scalar flelds in the dynamical 
evolution of matter and geometry has been studied self-consistently, mainly associated with the 
dynamics of non-equilibrium phase transitions!^, ^ in models that fall, broadly, in the 'small fleld' 
category. The conclusion of these studies is that a treatment of the quantum fluctuations that couple 
self-consistently to the dynamics of the metric provides a solid quantum fleld theoretical framework 
that justifies microscopically the picture based on classical inflation. At the same time these 
studies provide a deeper understanding of the quantum as well as classical aspects of inflation and 
inflationary perturbations. They clearly reveal the classicalization of initial quantum fluctuations!^, 
^, and furnish a microscopic explanation (and derivation) of the effective, homogeneous classical 
inflaton @. 

The purpose of this article is to provide a quantum treatment of models whose classical coun- 
terpart are large fleld models. The classical description in these models begins with a homogeneous 
inflaton scalar with very large amplitude (j) ~ Mp/||3|, |^, ^, i.e, very high up in the scalar potential 
well. The question of how is this initial condition achieved is typically answered in terms of a 
probabilistic distribution of initial conditions ^, as is the case for chaotic inflation. 

Instead, in this article we study the dynamics that results from the evolution of a quantum 
state which drives the dynamics of the scale factor through the expectation value of the energy 
momentum tensor. 

The initial state the we consider is a squeezed state with a large number of particles distributed 
in a shell in momentum. While squeezed states had been studied in quantum optics|lC] and also in 
cosmology [0] the state considered in this article is similar to that invoked to model the evolution 
of highly excited initial quantum states to describe the dynamics of heavy ion collisions. These 
states are characterized by a large population of quanta localized in momentum shells and had been 
coined tsunami-waves in[12, 13, 14|. The time evolution of these states leads to the formation of 
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a non-equilibrium plasma and features properties similar to those expected to occur in the cooling 



and expansion of a quark gluon plasma after the collision [ 15 1. 
The goals of this article: 



We here adapt the ideas and concepts in refs.|12, 13, 14] to study the self-consistent dynamics 



of the metric and the evolution of a highly excited quantum state with the goal of providing a 
quantum description of large field inflationary models without assuming an expectation value for 
the scalar field. We introduce novel quantum states, which are the cosmological counterpart of the 



tsunami-waves introduced in[Q3, 13, 14| with the following properties: 



Pure states: the states under consideration, defined by the wave-functional of the form 
are pure states. Other proposals involving a mixed state density matrix as initial state are 
discussed in sees. Ill, IV and the Appendix. 

Vanishing expectation value of the scalar field: Unlike most models of classical chaotic inflation 
in which the scalar field obtains an expectation value, taken as a classical field, the expectation 
value of the scalar field in the tsunami- wave states given by the wave- functional ( [l8| ) vanishes. 
[See sec. IV for a non-zero expectation value of the scalar field]. 



• Highly excited initial modes: the tsunami- wave state described by eq. (18) with the covariance 
kernel given by eqs. (|2^), (^) describes a state for which the modes inside a band are occupied 
with a non-perturbatively large 0(1/A) number of (adiabatic) quanta. This requires the use 
of non-perturbative methods as the large approach and makes the states considered to 
be very far from the vacuum. We remark that very high energy modes, those that will 
become super horizon during the last 10 or so e- folds, hence are of cosmological importance 
today, must be in the vacuum state so as not to lead to a large amplitude of scalar density 
perturbations ]P, ^, 

This type of quantum states is clearly a novel concept, it presents an alternative to typical infla- 
tionary scenarios that invoke the dynamics of a classical scalar field which in most cases ignore the 
quantum dynamics. 

The quantum nature of the tsunami states is due to the coherence between different modes 
and gives rise to dynamical consecuences. (See |14(]). In the present case the system is effectively 
classical only after the redshift has assembled the modes into a zero-mode effective condensate. 

The tsunami-wave states described above are the simplest states and will be the focus of our 
study. These states can be generalized to describe mixed-state density matrices and to also allow 
for an expectation value of the scalar field. These generalizations are described in sections III.C 
and IV and in the appendix and are found to lead qualitatively the same features revealed by the 
simpler pure states. 

We establish the conditions under which such quantum state leads to inflationary dynamics and 
study in detail the self-consistent evolution of this quantum state and the space-time metric. 

We emphasize that we are not proposing here yet a new model of inflation. Instead we focus 
on inflation driven by the evolution of a quantum state, within the framework of familiar models 
based on scalar flelds with typical quartic potentials. This is in contrast with the usual approach in 
which the dynamics is driven by the evolution of a homogeneous classical field of large amplitude. 

Brief summary: We find that inflation occurs under fairly general conditions that are the 
quantum equivalent of slow-roll. There are two consecutive but distinct inflationary stages: the 
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first one is completely determined by the quantum features of the state. Even when the expectation 
value of the scalar field vanishes at all times in this quantum state, the dynamics of the first stage 
gives rise to the emergence of an effective classical homogeneous condensate. The amplitude of the 
effective condensate is non-perturbatively large as a consequence of the non-perturbatively large 
number of quanta in the band of excited wavevectors. The second stage is similar to the familiar 
classical chaotic scenario, and can be interpreted as being driven by the dynamics of the effective 
homogeneous condensate. The band of excited quantum modes, if not superhorizon initially they 
cross the horizon during the first stage of inflation, hence they do not modify the power spectrum 
of scalar density perturbations on wavelengths that are of cosmological relevance today. Actually, 
in the explicit examples worked out here, the excited modes are initially superhorizon due to the 
generalized slow-roll condition. Therefore, in a very well defined manner, tsunami quantum states 
provide a quantum field theoretical justification, a microscopic basis, for chaotic inflation, explaining 
the classical dynamics of the homogeneous scalar field. 

In section II we introduce the quantum state, obtain the renormalized equations of motion for 
the self-consistent evolution of the quantum state and the scale factor. In section III we provide 
detailed analytic and numerical studies of the evolution and highlight the different inflationary 
stages. In section IV we discuss generalized scenarios. The summary of results is presented in the 
conclusions. We derive the equations of motion for mixed states in the appendix. 



2 Initial state and equations of motion 



As emphasized in the introduction, while most works on inflation treat the dynamics of the inflaton 
field at the classical level, we use a quantum description of the inflaton. 

We focus on the possibility of inflation through the dynamical quantum evolution of a highly 
excited initial state with large energy density. Consistently with inflation at a scale well below the 
Planck energy, we treat the inflaton fleld describing the matter as a quantum fleld whereas gravity 
is treated semiclassically. 

The dynamics of the classical space-time metric is determined by the Einstein equations with a 
source term given by the expectation value of the energy momentum tensor of the quantum inflaton 
fleld. The quantum field evolution is calculated in the resulting metric. 

Hence we solve self- consistently the coupled evolution equations for the classical metric and the 
quantum inflaton field. 

We assume that the universe is homogeneous, isotropic and spatially fiat, thus it is described 
by the metric, 

ds"^ = dt^ - a^{t) dx^ . (1) 

Anticipating the need for a non-perturbative treatment of the evolution of the quantum state, 
we consider an infiaton model with an A^-component scalar infiaton field <^(x) with quartic self- 
coupling. We then invoke the large N limit as a non-perturbative tool to study the dynamics]^, 
^, 12, 16 1 . This choice is not only motivated by the necessity of a consistent non-perturbative 



treatment but also because any grand unified field theory will contain a large number of scalar 
fields, thus justifying a large N limit on more physical grounds. 
The matter action and Lagrangian density are given by 



j d^X ^rn = j 



d^x a^{t) 



2 ^ ' 2 a^{t) 



V{^{x)) 



4 



m 



A 



F($) = — $^ + — ($^) +-^7^<I>^ 



1 



(3) 



and we will consider > 0, postponing the discussion of the case < to subsequent work. 
Here TZ{t) stands for the scalar curvature 

7^(^) = 6 



a{t) a?{t) 



(4) 



The ^-coupling of to the scalar curvature 'R-{t) has been included in the Lagrangian since it 

is necessary for the renormalizability of the theory. 

The discussion of the alternative inflationary mechanism that we are proposing and the physical 
description of the quantum states becomes more clear in conformal time 



dt' 
W) 



in terms of which the metric is conformal to that in Minkowski space-time 

ds"^ = a^{r) {dT'^ - dx^) . 
We introduce the conformally rescaled field 

*(r,x)=a(i)$(t,x) 
in terms of which the matter action becomes 



dTd^x < - [{Mf - (V*)^] - a^(r) V 



a{r) 



(5) 

(6) 

(7) 

(8) 



Since wc are interested in describing the time evolution of an initial quantum state, wc pass on to 
the Hamiltonian description in the Schrodinger representation. This procedure begins by obtaining 
the canonical momentum conjugate to the quantum field, n(T, x), and the Hamiltonian density 

n{r,-K) 



^'Cr,x), 



n(r,x) 
7^(r,x) 

H{T) = J d3xH(r,x) , 



a{T) 



(9) 



where the prime denotes derivative with respect to the conformal time T. 
In the Schrodinger representation the canonical momentum is given by 

S 



n»(r,x) 



(5*"(r,x) 



a = l,...,N . 



(10) 



The time evolution of the wave-functional T 
equation 



is obtained from the functional Schrodinger 



H 



d_ 



T 



(11) 
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The implementation of the large N limit begins by writing the field as follows 

^(x, T) = (ct(x, T), vf(x, T)) 

= (^/ivv'(r) + x(x,r),7f(x,r)), 



(12) 



where we choose the 'I'-axis in the direction of the expectation value of the field and we collectively 
denote by tt the — 1 perpendicular directions. 



(^(x, T)) 



(a(x, T)) 
(X(x, T)) = , 



(13) 



where the expectations value above are obtained in the state represented by the wave-functional 



T 



introduced above. 

The leading order in the large limit can be efficiently obtained by functional methods (see 
refs.|^, ^, |l^, ^ |16[ and references therein). The contributions of x to the equations of motion are 
subleading (of order 1/A^) in the large limit [P, 16 1. 

It is convenient to introduce the spatial Fourier modes of the quantum field 

7rfc(T) = I d^x7Ti^,T)e'''-^ (14) 
In leading order in the large limit, the explicit form of the Hamiltonian is given by|^, |9|, 14, |l^ 

H{T) = NVKiiT)-^ (Y^i^k-^-k)] +Y.H,iT), 

\ k / k 



Ki{r) 

Hk{T) 



i ^'2(T) + ^ ^2 V'2(r) + ^ v'(r), 



m + 4 7c H — 



7^(r) 



2 a'^{T) 2 Na?{T) 



(15) 

(16) 
(17) 



where V is the comoving volume. We assume spherically symmetric distributions in momentum 
space. 

That is, in the large A^ limit the Hamiltonian operator (^) becomes a time dependent c-number 
contribution plus a quantum mechanical contribution, Hk{T), given by a collection of harmonic 
oscillators with time-dependent frequencies, coupled only through the quantum fluctuations (vffc • 

In eqs.(|l5|)-([T7|) the scale factor a(T) is determined self-consistently by the Einstein-Friedmann 
equations. 
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2.1 Tsunami initial states 



To highlight the main aspects of the inflationary scenario proposed here, and to establish a clear 
difference with the conventional models, we now focus our discussion on the case of vanishing 
expectation value, i.e, ip{T) = 0, and a pure quantum state. The most general cases with mixed 
states described by density matrices and non-vanishing expectation value of the field are discussed 
in detail in sections IIIC and IV and in the appendix. 

For ip{T) = the quantum Hamiltonian (jl^) becomes a sum of harmonic oscillators with time 
dependent frequencies that depend on the quantum fluctuations. Therefore, we propose a Gaussian 
wave-functional of the form 



T 



1';T 



Afr{T) Yl . (18) 



The functional Schrodinger equation (|ll| ) in this case leads to evolution equations for the normal- 
ization factor A/t (T) and the covariance kernel (T) whose general form is found in the appendix 
(see also [^). The evolution of the normalization factor is determined by that of Ak, while the 
equation for Ak is 

iAi{T)=Al-u;l{T), (19) 

where primes refer to derivatives with respect to conformal time. As described in the appendix for 
the general case, the above equation can be linearized by defining (see appendix) 

where the mode functions ipk satisfy the equation 

^l+u;l{T)ipk = (21) 
In terms of these mode functions the self-consistent expectation value 

— [TTk ■ TT-k) 



N N J (27r)3 

(TTfc-Tr.fc) = — = — Iv'fel ■ (22) 

^AR,k ^ 

We now must provide initial conditions on the wave functional to completely specify the dynamics. 
Choosing the initial (conformal) time at T = with a(T = 0) = 1, the initial state is completely 
specified by furnishing the real and imaginary parts of the covariance A^ at the initial time. We 
parameterize these asQ 

AR,k{0) = nt ; Aj^kiO) = Wfe(O) 6k (23) 
Choosing the Wronskian of the mode functions 'Pk{T) and its complex conjugate to be 

"Pk ^k -'p'k'Pk = 2« (24) 



^In the case for which luI{0) < we choose ujk{0) = v'PT1A?2(0)~^;(0)76| 



7 



determines the following initial condition on the mode functions (see appendix) 



V9fc(0) = -= ; (/p'fc(O) = - [uJk{0)Sk + i ^k\ '^fc(O) . (25) 



where we also used eqs.(pO|) and (p3[). 

An important alternative interpretation of these mode functions is that they form a basis for 
expanding the Heisenberg field operators (solution of the Heisenberg equations of motion) 

^ [a, V'fc(r) ^"^-^ + 4 ^^(T) e-^-^] , (26) 
with (Xfc; o|. annihilation and creation operators, respectively, with canonical commutation relations. 



The Wronskian condition ( |24| ) ensures that the Tf{x,T) fields and their conjugate momenta obey 
the canonical commutation relations at equal conformal times. 

The physical interpretation of these initial states is highlighted by focusing on the occupation 
number of adiabatic states as well as on the probability distribution of field configurations. 

• Occupation number: it is at this point where the description in terms of conformal time proves 
to be valuable. In conformal time the Hamiltonian in the large A'^ limit is that of a collection 
of harmonic oscillators with time dependent frequencies. It is then convenient to introduce 
the adiabatic occupation number operator 



Mr) = ^ 



LOk{T) 2 



(27) 



with Hf: given by eq. (|T^). 

In particular the occupation number at the initial time is given by (see appendix) 

Here, the special case with 0^ = uJk{0) and 5^ = corresponds to the adiabatic vacuum 
(jifc = ). Instead, we study an initial state in which a band of wave- vectors are populated 
with a non-perturbatively large number of particles. More precisely, we consider initial states 
for which 

n. = o(i) inside .he excited band, 

TT-fc = outside the excited band. . (29) 

where A is the quartic self-coupling. 
This is accomplished by choosing, 

O ( TT I inside the excited band. 



— = — -— and (5fc = outside the excited band. (30) 
SZfc cjfc(O) 
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These initial states are highly excited, the expectation value of the energy momentum tensor 
in these states leads to an energy density ~ 1/A and are, therefore, non-perturbative. We 
will refer to the case where the excited band is narrow as the narrow tsunami. 

It must be stressed that the particle distribution alone partially determines the initial 
state. As we see from eq.(^), the initial state is completely defined specifying two functions 
of A; : rifc and 5k ■ 

Probability distribution: an alternative interpretation of these initial states is obtained by 
focusing on the probability distribution of field configurations at the initial time. It is given 
by 

V[n] = T T = oj ^ = AA(0) ]J e"^* ^^'^-'^ (31) 

k 

An intuitive quantum mechanical picture of the wave-functional for the modes in the excited 
band is the following. At the initial time the instantaneous Hamiltonian corresponds to a 
set of harmonic oscillators of frequencies ujk{0), while the width in field space of the initial 
Gaussian state is determined by $7^ . For a mode in the vacuum state l/\/TTk ~ l/y'uJkiO) 
and the typical amplitudes of the field are vTfc ~ 1/ y^u!k{0) which is the typical width of 
the potential well. While for a mode inside the excited band the width in field space 
is ~ LL'fc(O) [see eq.(|30|)]. Thus, large amplitude field configurations with 

T^k^ko ~ 1/y^A ajfc(O) » l/y^ujfjfi) have a probability of 0(1), i.e, large amplitude config- 
urations within the band of excited wave- vectors are not suppressed. That is, the width of 
the probability distribution for these modes is much larger than the typical size of the poten- 
tial well and there is a non-negligible probability for finding field configurations with large 
amplitudes of 0(1/ A). 

These highly excited initial states had been previously proposed as models to describe the 



initial stages of ultrarelativistic heavy ion collisions and had been coined 'tsunami waves' [12 
p^ , p^ ] . They represent spherical shells (in momentum space) with large occupation numbers 
of quanta, describing a state with a large energy density with particles of a given momentum. 



2.2 Back to comoving time: renormalized equations of motion 



Having set up the initial value problem in terms of the tsunami- wave initial wave-functionals, the 
dynamics is now completely determined by the set of mode equations eqs.(21) with eqs.(^6|)-(^) 
and the initial conditions eqs.(^5|). However, in order to establish contact with more familiar results 
in the literature, it is convenient to re-write the equations of motion in comoving time. This is 
accomplished by the field rescaling given by eq.(^) which at the level of mode function results in 
introducing the comoving time mode functions fkit) related to the conformal time ones '^k{T) as 



fk{t) 



a{t) 



The equations of motion in comoving time for these mode functions are 



fk{t)+3H{t)fkit) + 



a^t) 



M\t) 



fk{t) = 



(32) 



(33) 
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m' + mt) + ^ I (^l/^WP (34) 

/fc(0) = ^; fk{0) = -[uJk{0)Sk + H{0)+ink]fk{0). (35) 
V"fe 



The Einstein-Friedmann equation are, 

«^<"Kly)^^^ • '''' 

where the expectation value is taken in the time evolved quantum state. It is straightforward to 
see that the expectation value of the energy-momentum tensor has the perfect fluid form, as a 
consequence of the homogeneity and isotropy of the system |8|, |9[. 

Thus the set of equations (^)-(|3^) provide a closed set of self-consistent equation for the 
dynamics of the quantum state and the space-time metric. 

2.2.1 Renormalized Equations of Motion in the Large N limit 

The set of equations that determine the dynamics of the quantum state and the scale factor need 
to be renormalized. The field quantum fluctuations 



N J 2(27r) 

requires subtractions which are absorbed in a renormalization of the mass, coupling to the Ricci 
scalar and coupling constant. The expectation value of the stress tensor also requires subtractions 
(but not multiplicative renormalization). Since the divergence structure is determined by the large 
energy, short distance behavior, the band of excited modes does not influence the renormalization 
aspects. Therefore, we use the extensive work on the renormalization program which is available 
in the literature referring the reader to references|^, |9| for details. We here summarize the aspects 
that are most relevant for the present discussion. 

First, it is convenient to introduce the following dimensionless quantities, 

,/ N H(t) k 

T = m t ; n[T) = ; q = — ; 

m m 



- ; % = — ; 9 = ^ ; /.(r) = A(t) , (37) 



where m and A stand for the renormalized mass of the inflaton and the renormalized self-coupling, 
respectively In terms of these dimensionless quantities we now introduce the dimensionless and 
fully renormalized expectation value of the self-consistent field as 



(38) 



where the terms subtracted inside the integrand renormalize the mass, coupling to the Ricci scalar 
and the coupling constant & M. The dimensionless and renormalized expressions for the energy 
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density e and pressure p are given by 
A 



6(t) 



2N m4 



/T^00\ 



+ 



+ ? / 9' dq\\f,{rr-S^{q,T) + 



p{t) 



A 



2iV m4 



4 2 



a2(r) 



[|/,(r)|2-S2(g,T)] 



[\Mr)\'-S2{q,T)] 



(39) 



Where the renormahzation subtractions 5i and S2 are given by,| 



S2{q,T) 
B(t) 



q 



+ 



1 



a4(r) 2ga4(r) 



[B{t) + 2d- 



a21 



+ 



e(g-i) 

8g3 a4(r) 
1 



-S(r)2 - a{TfB{T) + 3a(T)d(r)S(r) - 'ih^{T)B{T) 



a2(r)[l+5S(r)] . 



B(t) + 



Q(g-i) 

8g5 a2(r) 



3S(r)2 + a(r)-^ \a{T)B{T) 

ClT I 



(40) 



We choose here = (minimal couphng), the renormahzation point k = \m\ and a(0) = 1. 

In summary, the set of coupled, self-consistent equations of motion for the quantum state and 
the scale factor are 

„2 



+ l + (7S(r) 



U{r) 







/,(0) 



1 



; /q(0) = -K(5, + /i(0)+iJ^,]/,(0) 



lq^ + 



1 + 5^(0) 



no) 



6m? 



plus the Einstein-Friedmann equation of motion for the scale factor 



hHr)=L^eir) 



where = 



3M|,;A 



(41) 
(42) 

(43) 
(44) 



with 5S(t) and e(T) given by eqs. ( pq) and (39) respectively. 

In order to implement the numerical analysis of the set of eqs. (pll)-(^), (]38| ) and (^) we 
introduce an ultraviolet momentum cutoff A. For the cases considered in this article we choose A ~ 
200 and found almost no dependence on the cutoff for larger values. As befits a scalar inflationary 
model, the scalar self-coupling is constrained by the amplitude of scalar density perturbations to 
be A ~ 10~^2[^, ^] implying that g < 10~^^. Therefore, the subtractions can be neglected because 
S, ~ 0{gA^) < 10-4. 

The initial state is defined by specifying the 0,q and 6q. We determine the range of these 
parameters and 5q by the excitation spectrum for the tsunami-wave initial state, as well as the 
condition that lead to infiationary stage. This will be studied in detail in the next section. 
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3 Tsunami inflation 



As emphasized in the previous section, the scenario under consideration is very different from the 
popular treatments of inflation based on the evolution of classical scalar inflaton field ||2|, ^, ^, 
In these scenarios all of the initial energy is assumed to be in a zero mode (or order parameter) at 
the beginning of inflation and the quantum fluctuations are taken to be perturbatively small with 
a negligible contribution to the energy density and the evolution of the scale factor. 

In contrast to this description, our proposal highlights the dynamics of the quantum states as 
the driving mechanism for inflation. The initial quantum states under consideration correspond to 
a band of quantum modes in highly excited states, thus the name 'tsunami-wave'[12, 14 1. This 



initial state models a cosmological initial condition in which the energy density is non-perturbatively 
large, but concentrated in the quanta rather than in a zero mode. 

We now study under which general conditions such a state can lead to a period of inflation 
that satisfies the cosmological constraints for solving the horizon and entropy problems entailing 
the necessity for about 60 e-folds of inflation. 

It is understood that inflation takes place whenever the expansion of the universe accelerates, 

i.e, 

- = /i2 + A = - — [e + 3p] >0, (45) 

with L given in eq. (^^ and e and p given by eqs. (^9]). 

While our full analysis rely on the numerical integration of the above set of equations, much we 
learn by considering the narrow tsunami case. 

3.1 Analytical study: the narrow tsunami case 

Before proceeding to a full numerical study of the equations of motion, we want to obtain an 
analytic estimate of the conditions under which a tsunami initial quantum state would lead to 
inflation. 

Our main criterion for such initial state to represent high energy excitations is that the number 
of quanta in the band of excited modes is of 0{1/X). This criterion, as explained above, is tanta- 
mount to requiring that fleld configurations with non-perturbative amplitudes have non-negligible 
functional probability. Progress can be made analytically by focusing on the case in which the band 
of excited field modes is narrow i.e, its width A/c is such that A/c ^ or in terms of dimensionless 
quantities Aq/qo <C 1. We introduce the following smooth distribution 

n, = -J , with ^«l, (46) 

g 

with ujq given by eq. (^3|) and A/q a normalization constant that fixes the value of the total energy. 
In addition, we choose 6q = —h{0)/u!q as we discuss below in eq. (^7|) . 

This initial distribution posses the main features of the tsunami state described in the previous 
section. Since 51 <C 1, we have for q ^ qo, 

7^~->l ~->l (47) 

9 9 
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corresponding to highly excited states. While for \q - 



> Aq 



1 1 

— ~ — =^ ng ~ . 



(48) 



Thus, these modes are in a quantum state near the conformal (adiabatic) vacuum at the initial time, 
with Uq the number of quanta defined by eq. (^8|) in terms of dimensionless variables. For these 
distributions (narrow tsunamis), the integral over mode functions for the quantum fluctuations 
5(S(t) [given by eq. (|38|)] is dominated by the narrow band of excited states with mode amplitudes 
~ 1/ and can be approximated by 



5S(t) = g Aq ql |/,„(t)|2 + 0{g) + 0{g Aq) ^ |</.,„(t)|2 , 
where we have introduced the effective qo mode 

4>qo{T) = Va Aq qo fqoir). 



(49) 



(50) 



we note that the initial condition ( |4^ ) and the tsunami- wave condition (47) entail that despite the 
presence of the coupling constant in its definition, the amplitude of the effective qo mode is of 0(1). 
The equation of motion for the effective g^o-^ode takes the form 



(?!>g„(r) +3/i(r)(/.y„(T) + 
The scale factor follows from 



h\T)=L^ e(r 



+ l + l</'go(T)| 



(51) 
(52) 



with energy and pressure. 



{p + e){T) 



(r)P + 



3a2(r) 



(r)P + 
I<A,o(t)| 



^90 



{r)\' + 



2a2(r) 



(53) 



where we have neglected terms of 0{g). 

We will refer to the set of evolution equations (|5l|)-(|53[) as the one mode approximation evolution 
equations. 

In particular, within this one-mode approximation, the acceleration of the scale factor obeys 



a(r) 
a(r) 



^qa\ 



Therefore, the condition for an inflationary epoch, a > 0, becomes 



'^qoy 



< 



1 



2 ■ 4 

A sufficient criterion that guarantees inflation is the tsunami slow roll condition 

\'Pqoi^)\ « I 090 W I 



1 



'^qQ\ 



(54) 



(55) 



(56) 
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The initial conditions (^2|) and the condition that ~ 5 ^ 1 imply that the tsunami slow roll 
condition 



at early times is guaranteed if 5„o is such that 



(57) 



Hence tsunami- wave initial states that satisfy the tsunami slow-roll condition (57) lead to an infla- 
tionary stage. 

Moreover, in order to have slow roll (^) at later times, the effective friction coefficient 3/i(t) 
should be larger than the square of the frequency in the evolution equations (^). That is. 



gg + l + ffS(0) 
3/i(0) 



<. 1 . 



(58) 

implies that 



(i.e. the qo-mode should be deep inside the overdamped oscillatory regime). Eq.(^ 
/i(0) S> 1 and this together with eq.(|57|) implies that 6qg must be negative. 

A remarkable aspect of the narrow tsunami state is that it leads to a dynamical evolution 
of the metric similar to that obtained in classical chaotic inflationary scenarios in the slow roll 
approximation In particular the expression for the acceleration ( |5^ and the tsunami slow 

roll condition (|5€) are indeed similar to those obtained in classical chaotic inflationary models 
driven by a homogeneous classical field (zero mode). However, despite the striking similarity with 
classical chaotic models, we haste to add that both the conditions that define a tsunami state 
and the tsunami slow roll condition ( |5^ ) guaranteed by the initial value ( |57[ ) is of purely quantum 
mechanical origin in contrast with the classical chaotic-slow-roll scenario. Furthermore, we recall 
that the expectation value of the scalar field vanishes in this state. 



3.1.1 Early time dynamics 



Under the assumption of a tsunami wave initial state and the tsunami slow-roll condition (56) the 
contribution 0go(r) in the energy and in the pressure [see eq. (|53|)] can be neglected provided. 



3a2(r) 



\<P<ioir)\' 



(59) 



We call TA the time scale at which this rely no longer holds. Furthermore, we can approximate 
090 (^) by (I)goi0) if 

' (0) 



TA < 



(0) 



(60) 



This condition is fulfilled at least for ^ 1 due to the tsunami slow-roll condition ( |56| 
During this interval the Friedmann equation (^2|) takes the form. 



\a{r)] 


2 









D 



a2(r) 



+ E . 



(61) 



where we used that gS(0) = |i;^>ijg(0)p. This equation is valid as long as the characteristic time scale 
of variation of the metric is shorter than that of the mode (pq^ (r) . 
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The preceding equation can be integrated with solution 



sinh ( \fE T + c 



a{T) = 
Kr) = 

where the constants D, E and c are given by, 



a(r) 

E ^ ■ " ■ > 

coth (V^r + c) , A(r) = 



E > Q 



E 



sinh^ ( \/E t + c 



(62) 



sinh ~ y ^ ' 



(63) 



and 5So = 5S(0) = |</>,o(0)|2. 

We see from eqs. ( |6^ ) that during this interval there is an inflationary stage with an accelerated 
expansion = E > 0. We also see that /i(t) decreases with time until it reaches the constant 

value ^/E that determines the onset of a quasi-De Sitter inflationary stage. 

We now estimate the range of validity of the solution in eqs. (]6^). The flrst condition in eq. ( ^9|) 
is more stringent than the second one in eq. (|60|) for most of the interesting range of parameters. 



TA determines the time scale at which the solution (62) ceases to be valid, and the condition 
|59|) is no longer fulfllled, i.e. 



I^o(^a)P 



3a2(r^) 



(64) 



When this equation is valid, we see from eq. ( |53D that p + e becomes of the order of |(/)gQ(ryi)p. 
Furthermore, the above condition together with the slow roll condition, eq. (^), leads to 



< 1 . 



(65) 



Therefore, from eq.(|6l|) we see that /i(r^) ~ vE and is slowly varying. 

Since the slow roll condition guarantees that |(?!'g(T^)| <C |(/>g(ryi)|, we can now use eq.(^9|) along 
with the evolution equation (^) which setting h = \fE leads to the following relation 



From eq.(|6^) the scale ta is therefore estimated to be given by 



(66) 



TA 



ArgSinh 



VSqoE 





1 


)-] 






~ Ve 



(67) 

This initial inflationary period with a decreasing Hubble parameter exists provided the r.h.s. is 
here positive, i.e., 

,0 > . (68) 



/3E 
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In order to distinguish this phase from the later stages, to be described below, we refer to this early 
time inflationary stage as 'tsunami-wave inflation' because the distinct evolution of the scale factor 
during this stage is consequence of the tsunami- wave properties. 
A.i T = TA we have: 

^irA)^^^^ , h{rA)^VE. (69) 

For T > TA, qo/a{T) « 1 and the physical wavevectors in the excited band have red-shifted 
so much that all terms containing qq become negligible in the evolution equations. Therefore, all 
modes in the excited band evolve as an effective g = mode. Hence for t > ta the dynamics of the 
scale factor is described by an effective homogeneous zero mode and describes a different regime 
from the one studied above. Such regime is akin to the classical chaotic scenario. 



3.1.2 The effective classical chaotic inflationary epoch 

For T >> TA, when q^/a^ <C <C 1 all the physical momenta corresponding to the 

comoving wavevectors in the excited band have redshifted to become negligible in the equations of 
motion. The dynamics is now determined by the following set of equations for the effective zero 
mode and the scale factor, 

<^,„ (r) + 3 /i(r) (r ) + [1 + I (r ) 1 2] (r) = , 

h\T) = e(r) , (70) 
where the energy and pressure are given by, 



<r) = l|0,„(r)p + i|0,„(r)|2 + i|</>,„(r)r 



{p + e){T) = \^,Ar)?- (71) 

The initial conditions on (pq^^ and 0q„ are determined by their values at the time ta, while the 
slow-roll condition (^) determines that the imaginary parts of and (^^q are negligible. 

Therefore, after ta the dynamic is identical to that of a classical homogeneous field (zero mode) 

r/e//(r)=Re[(A,o(r)], (72) 

that satisfies the equations of motion, 

^e// + 3 /ir)e// + (1 + ^/e//) = > 

h\T) = L'e{T). (73) 

with energy and pressure. 



e r 



1-2 2 

2 ^e// + 2 ^''ff 4 ^^^^ 
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(p + e)(T) = vlfj, (74) 
and initial conditions [using eq. (|69|)], 



VeffiTA) = <t)qo{TA) = (/'qo(O) = ' 

VeffiTA) = Ki'^A) = - 3^(^^° "leffiTA) ■ (75) 

Where the value of rjeffiTA) is determined by the slow roll condition (^ </)g(,(ryi) ~ 0), the evolution 
eq. (^) and eq. (|65| ) and a(r^) and /i(ta) are given by eq. (|69|) . 

When (^Sq ^ 1, the quadratic term in the potential dominates, and we can integrate the 
previous equations to obtain 

^e//(r) = 7?e//(rA) - — (r - ta) , (for ^Sq < 1) . (76) 

This evolution is similar to that of classical chaotic inflationary models 1^, ^. Therefore for r > 
TA when the physical momenta in the excited band have redshifted so much that their contribution 
in the equations of motion of the quantum modes and the energy and pressure become negligible, the 
evolution of the quantum modes and the metric is akin to a classical chaotic inflationary scenario 
driven by a homogeneous c-number scalar field. This equivalence allows us to use the results 
obtained for classical chaotic inflation. Thus, as the classical slow roll condition (|?7e//| ^ \Veff\) 
holds, the evolution of the effective scalar field is overdamped and the system enters a quasi- 
De Sitter inflationary epoch. This inflationary period ends when the slowly decreasing Hubble 
parameter becomes of the order of the inflaton mass, i.e, 3 /i ~ 1 + i]eff ^^^^ stage the effective 
classical field exits the overdamped regime and starts to oscillate, the slow roll condition no longer 
holds and a matter dominated epoch {\rieff \ ~ l^e//| =^ p ~ 0) follows. 

However, we emphasize that while the effective zero mode ?7e// obeys a classical equation of 
motion and that the components of the energy momentum tensor eq.(|7^ are those from a classical 
field, the origin of this mode is purely quantum mechanical. From the identification (|75| ) and eq. 
(|4^ it is clear that the effective zero mode is a collective superposition of modes in the highly 
excited band. From the initial and tsunami wave conditions fq^irA) ~ l/y^H^ ~ it follows 

that the amplitude of the effective zero mode is r]f.ff{TA) ~ qo\/Aq. Restoring the dimensions and 
the proper powers of the coupling that were absorbed in the constant L in the Friedmann equation 
we find that the equations of motion ( [73|) with the stress tensor components (|7^ ) are those obtained 
from a (dimensionfull) classical homogeneous field ipeff{t) with a classical potential V{ipeff) given 
by 

feffit) = Veff{T) 

^i^eff) = -Y ^Iff + 8 V^f// (77) 
with the initial value at the time Ia = TA/m 
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The non-perturbative amplitude of the effective zero mode is a consequence of the non-perturbative 
amphtude of the excited quantum modes with an 0{1/X) number of quanta. 



3.1.3 Number of e-folds 

An important cosmological quantity is the total number of e-folds during inflation. As discussed 
above, there are two different inflationary stages, the first one is determined by the equations 
(|6l|)-(|62|) and characterized by a rapid fall-off of the Hubble parameter approaching a quasi-De 
Sitter stage. This new stage has been referred to as the tsunami-wave inflationary stage above 
to emphasize that the dynamics is determined by the distinct characteristics of the tsunami-wave 
initial stage. 

The second stage is described by an effective zero mode and the evolution equations 
and is akin to the chaotic inflationary stage driven by a classical homogeneous scalar field. The 
crossover between the two regimes is determined by the time scale ta (in units of the inflaton mass) 
and given by eqs.(67) at which the contribution from the term q^/a?'{T) to the equations of motion 



becomes negligible. Therefore there are to distinct contributions to the total number of e-folds, 
which is given by 

iVe('7o, /i(0)) = log a{TA) + iVe(0, Kta)) , (79) 

where a{TA) is given by eq. (^) and A'^e(0, h{TA)) is just the number of e-folds for classical chaotic 
inflation with an initial Hubble parameter /i(r^). 

We can express h{TA) as a function of qo and h{0), 



Hta) = M/f^ + ^ = V^^(O) - f 9^0 . (80) 



The number of e-folds during the first stage, is given by 



loga(r_A) ~ log 



1 + 5^0 



(81) 



The expression for the number of e-folds during the following, chaotic inflationary stage simplifies 
when (jrSo ^1. In this case the quadratic term in the potential dominates, and we can obtain 
simple analytical expressions 



iVe(0, hiTA)) = — r^iffirA) = — 5S0 = — , (for gl^o « 1) • (82) 



We see that the number of efolds grow when qo decreases at fixed initial energy eo- That is, we 
have more efolds when the energy is concentrated at low momenta. 



3.1.4 In summary 

Before proceeding to a full numerical study of the evolution we summarize the main features of the 
dynamics gleaned from the narrow tsunami case to compare with the numerical results. 

• The conditions for tsunami- wave inflation are i) a band of excited states centered at a 
momentum with a non-perturbatively large 0{l/g) number of quanta in this band, and 
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ii) the tsunami slow-roll condition eq.([5q). These conditions are guaranteed by the initial 
conditions on the mode functions given by eq.(E3) with the tsunami- wave distributions of the 



general form given by eqs. (Hq), (57). 



There are two successive inflationary periods. During the first one, described in sec. p.l.l] , 
the dynamics is completely characterized by the distinct features of the tsunami-wave initial 
state, the Hubble parameter falls off fast and reaches an approximately constant value ^/E 
that characterizes the quasi-De Sitter epoch of inflation of the second period. The second 



stage, described in sec. 3.1.2 can be described in terms of an effective classical zero mode 
and the evolution of this effective mode and that of the Hubble parameter are akin to the 
standard chaotic inflationary scenario. 

• The tsunami-wave initial state can be interpreted as a microscopic justification of the classical 
chaotic scenario described by an effective classical zero mode of large amplitude. The ampli- 
tude of this effective zero mode is non-perturbative as a consequence of the non-perturbative 
C(l/A) number of quanta in the narrow band of excited modes. Thus the initial value of 
the effective, classical zero mode that describes the second, chaotic inflationary stage, is 
completely determined by the quantum initial state. 

• An important point from the perspective of structure formation is that the band of excited 
wavevectors centered at qo either correspond to superhorizon modes initially, or all of the 
excited modes cross the horizon during the first stage of inflation, i.e, during the tsunami 
stage. This is important because the chaotic second stage of inflation which dominates during 
a longer period guarantees that the band of excited modes have become superhorizon well 
before the last 10 e-folds of inflation and hence cannot affect the power spectrum of the 
temperature anisotropics in the CMB. The fact that the tsunami-wave initial state is such 
that the very high energy modes (necessarily trans-Planckian) that cross the horizon during 
the last 10 e-folds and are therefore of cosmological importance today are in their (conformal) 
vacuum state leads to the usual results from chaotic inflation for the power spectrum of scalar 
density perturbations. 

Although these conclusions are based on the narrow tsunami case, we will see below that a full 
numerical integration of the self-consistent set of equations of motion confirms this picture. 

In sections 3.2 and Q we show how this results can be easily extended to more general particle 



distributions and more general initial states. 
3.2 Numerical example 

To make contact with familiar models of inflation with an inflaton field with a mass near the grand 
unification scale, we choose the following values of the parameters: 

10-^ , A = 10-12 , AT = 20 (83) 



Mpi 

where the number of scalar fields = 20 has been chosen as a generic representative of a grand 
unified quantum field theory. For these values we find 

= ^— = 3.35 • 10^ . 

3M|,;A 
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As an example we shall consider an initial energy density pp = (Too) = 10'^ Ml.i. Thus, 
the initial value for the Hubble parameter is Hq = ^/STrpo/SMpi = 3.53 • 10^^ GeV{= 1.654 • 
10^^ km/s/Mpc). These initial conditions in dimensionless variables give eg = 2.50 and h{0) = 
2890. 

In addition, the slow roll conditions (|58| ) imply: 



which in this case results in 



3/i(0) 



go < 95 . 



We choose qq = 80.0, and initial conditions in eq. (^) with Qq and 6q given by eq. (|4q ) and (pT\j. 
These initial conditions satisfy the tsunami slow roll condition, 

|wg <5q + /i(0)| « 1 (84) 

Furthermore, we take Aq = 0.1 and Mn is adjusted by fixing the value ^^(O) = gTiQ which for the 
values chosen for eo and by eqs.(p9|) and (^) and ( p6| ) gives ^Sq = 7.81 • 10"*^. 

Figure |^ displays eo vs go along lines of constant number of e-folds , while figures display 
the solution of the full set of equations (|4l|)-(44) with (|39|). An important feature that emerges 



from these figures is that for the set of parameters that are typical for inflationary scenarios and 
for large values of go = ko/m (but well below the Planck scale) the number of e-folds obtained is 
more than sufficient as shown by fig.^. 

We also show that the dynamics of the full set of equations (|4l|)-(44) with ( |39| ) is correctly 



approximated by the narrow tsunami case studied in the previous subsections: the one mode 
approximation [eqs. (|49|)-(|53|)], the early time analytical formulae (for r < ta) [eqs. (|62|)-(|6^)], and 
the effective classical field (for r > ta) [eqs. (0)]. The agreement between the analytic treatment 
and the full numerical evolution is displayed in figures |2[|^. 

The early time analytic expressions predict an inflationary period during which the Hubble 
parameter falls off fairly fast, that lasts up to ta ~ 0.133 [eq. (|67|)] reaching an asymptotic value 
of /i(ta) = 36.2 [eqs. (|63|) and (39)]. The one mode approximation gives the same prediction 



h{TA) = 36.1, and numerically evolving the full set of equations we flnd h{TA) = 35.5. Thus, we 
see from this values and from flgs. ^ and that both approximations are fairly accurate for early 
times. 

After TA, the geometry reaches a quasi-De Sitter epoch. We have shown in the previous sub- 
section that after the time ta the evolution equations for the one mode approximation reduce to 
those of an effective classical fleld. The effective zero mode approximation correctly predicts the 
dynamics in this epoch as can be gleaned from flgures |^|^. 

While the stage of early tsunami inflation up to ta results in only 8.5 efolds, the following 
quasi-de Sitter stage described by the effective classical scalar field lasts for a total of 1900 efolds. 
For the values of parameters chosen above, ^Sq <C 1, hence we can estimate the number of efolds 



with eq. (82). Using eq. ([79|) we obtain a total of 1970 efolds while the one mode approximation 



yields 1960 efolds. Both results agree with the full numerical solution of the equations (see fig. |6D. 

Furthermore, as stated above infiation ends when h ~ — 3^ ~ 3, after which a matter domi- 
nated epoch follows. 
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3.3 Other distributions and other states: 



The vahdity of the physical picture that emerges from the previous analytic and numerical study 
is not restricted to pure states or narrow distributions of the form given by (^). We have also 
studied more general distributions and mixed states: 

Other distributions: 

The narrow tsunami case where a single quantum mode qq dominates the dynamics has been 
extremely useful to study the dynamics in the previous section. The generalization to the case with 
continuous distributions of g-modes can be easily obtained making the changes: 



l'A,o(^)l' 



9 j q^dq\f,{T)W 
dq\U{T)\\ 
dq g'|/,(T)p. 



a / 9 



g 



(85) 



The two stages of inflation are always present for such continuous modes distribution as long as 
the following generalized slow-roll conditions is fulfilled 

(86) 



9 j q' dq\fgir)\' <.g J |/,(t)|^ 

that imposes on 5q the condition \ujq 5q + h{0)\ <C 1. 

The effective zero mode in the second stage of inflation is now given by 



9 J q^dq |/,(r) 



Our numerical study with general distributions reveals that the analytical picture obtained by 
substituting eq.([85|) in sec. IIIA correctly reproduce the dynamics. 

Other (mixed) states: Although we have focused for simplicity on tsunami pure initial states, 
we have also investigated the possibility of mixed states. Mixed state density matrices and their 
time evolution are discussed in the appendix. The mixing can be parametrized in terms of angles 
@k as given in equation (105) and the number of (conformal) quanta are given by eq. ( |119| ). The 
only relevant changes that occur are in the integrals for S(t) ; e(r) ; p{t) in which 



l/.Wl'-|/,(r)pcoth 



l/.(T)|'-|/,(r)pcoth 



(87) 



Tsunami and slow-roll conditions on the mode functions given by eqs. (|30|) and ( |57D lead to 
tsunami-wave inflation followed by chaotic inflation just as discussed above. In the narrow tsunami 
case the only change is that the effective go-mode is rescaled by the mixing factor, i.e, 



e 



<?0 



(r)|'-|0,o(r)|'coth 



e 



90 



It is also illuminating to contrast the tsunami-wave mixed states with the more familiar thermal 
mixed states. The latter are obtained by the choice 

T 



LOq ; 6q 



; e„ 



(89) 
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with T some value of temperature. In this case it is straightforward to see that (quantum) equipar- 
tition results in that the contributions of the modes and their time derivatives to the energy and 
pressure are of the same order (|/g(T)p ~ [/i(0)^ + ^g] l/g(''")P)- Hence, for these thermal mixed 
states the tsunami slow-roll condition is not fulfilled. This is obviously not surprising, such a choice 
of thermal mixed state leads to a FRW epoch which is not inflationary. Hence the tsunami-wave 
initial conditions along with the generalized slow-roll conditions lead to two successive inflationary 
epochs in striking contrast to the familiar mixed thermal states. 



4 Generalized chaotic inflation 

The previous analysis, confirmed by the numerical evolution of the full self-consistent set of equa- 
tions leads to one of the important conclusions of this article, that tsunami-wave initial states 
provide a microscopic justification of the chaotic inflationary scenario. 

We have focused our discussion on initial states with vanishing expectation value of the scalar 
field (order parameter) and where the energy is concentrated in a momentum band (tsunami initial 
states). This choice brings to the fore the striking contrast between this novel quantum state and 
the usual classical approach to chaotic inflation. In this section we study the dynamics in the 
case in which the initial state allows for a non-vanishing expectation value of the scalar field along 
with some of the initial energy localized in excited quanta. We refer to this case as generalized 
chaotic inflation to distinguish from the tsunami-wave state studied above. This generalization thus 
includes both cases: the classical chaotic inflation in the limit when there are no excited modes, as 
well as the tsunami initial state when all of the energy is localized in a band of excited modes and 
the expectation value of the field vanishes. 

The relevant equations of motion in comoving time for the mode functions in this case are 



given by (121) in the appendix. Along with the dimensionless variables (^) it is also convenient 



to introduce a dimensionless expectation value as 

ri'ir) = ^ <P\t) (90) 

In this generalized case with rj ^ the equations of motion for the mode functions /<j(t) (in terms 

of dimensionless variables) are the same as in eq. (^) after the replacement (7S(r) 5'S(r) +r/^(r) 
and the equation of motion for r/(r) is given by 

^ + 3Mr)^ + [1 + r^\r) + g^r)] r/(r) = (91) 

V{0)=V0 ; 7/(0) = r?o (92) 



The Einstein-Friedmann equation is given by ( [44D but with the energy density and pressure now 
given by 

+ |y"</'dg||/,|'-5i((/,T) + ^ [|/,|2-S2(q,t)]| , (93) 



(p + e)(T) = f +5y"9'd(/||/,|2-5i(g,r) + ^[|//-52((/,r)]| 



(94) 
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where the renormalization subtractions 5i ; 5*2 are obtained from those given by eq. ( |40| ) upon the 
replacement gT, — > gT, + rj"^. 

From this expression we see that for fixed (large) energy density Status: RO 
there are two different possibilities: if the zero mode squared r/^(0) is larger than the quantum 
fluctuations gT, and g J dq the dynamics is basically similar to that in the usual chaotic 

inflationary scenarios. This corresponds to most of the initial energy density to be in the zero mode 
and little energy density in the band of excited states. On the other hand, for small ?/^(0) most of 
the initial energy density is in the tsunami quantum state and the initial dynamics is akin to the 
r] = case. To quantify this statement and clarify the interplay and crossover of behaviors between 
the r] = and the generalized chaotic case, we now resort again to the narrow tsunami case, which 
highlights the essential physics. 



The relevant equations are: i) the equations of motion for the effective go~™ode (50), 

«2 



+ l+r,'{r) + \^,,{T 



and for the zero mode r] 

iiir) + 3/i(r)??(r) + [l + 7^\t) + |0,o(^)l'] v{r) = 
and the Hubble parameter given by (^2|) with the energy density given by 



e(T) = - 



{rW + v\r)]+- (|0,„(r)p + 7y2(^))+ (|<^^^(^)|2 + ^2(^))^ ^ 



2a2(r) 



(95) 



(96) 



(97) 



The acceleration of the scale factor is now given by 



a(r) 
a(r) 



koir)\' + f (r) - i {\^,,{T)f + rj^r)) - \ (|0,„(r)p + r^^r))' 



(98) 



where again we have neglected the renormalization contributions and terms of 0{g) consistently in 
the weak coupling limit g <^\. From eq. (p^) the generalized condition for an inflationary epoch 
(within the narrow tsunami case) becomes 



(r)p+f(r)< 



which is fulfilled if the following generalized slow roll condition holds 

■ ,(^)l'+^'W«l</'.o(^)l' + ^'W 



(99) 



(100) 



Under these conditions and from eqs.(p5|)-(p6|) and the dynamics of the scale factor driven by the 
energy density eq.(|97|) we can now distinguish the following different inflationary scenarios. 

• Tsunami dominated: When {q^ + 1) |</'gQ(0)p > '7^(0) the excited states in the tsunami- wave 
carry most of the initial energy density. In this case the results of the previous section apply 
and the scale factor takes the form as in eq.(^) with D given by eq.(|6^), and E given by 
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eq.(|63|) but with gTjQ gTjQ + ??^(0). There are two consecutive inflationary stages as in the 



previous section. The flrst described by eq.(|6^), lasts up to the time scale ta defined by 



at which the red-shift of the momentum qq is such that qo/a^TA) ^ 1- The secondary stage 
is a usual classical chaotic inflationary epoch determined by the dynamics of an effective zero 
mode given by 

Veff{r) = v\r) + \<P,,{T)\^ (101) 

because for t > ta, qQ/a^ir) <^ 1 and the effective equation of motion for (pq^ir) is the same 
as that for ??(t). 

• Zero mode dominated: When ?7^(0) ^ (qq + 1) |i;^go(0)p the energy density stored in the zero 
mode is much larger than that contributed by the excited states in the tsunami-wave. In 
this case the energy density eq.(^) is completely dominated by the zero mode. The ensuing 
dynamics is the familiar classical chaotic scenario driven by a classical zero mode, without 
an early stage in which the scale factor is given by eq.(|62[) which is the hallmark of the 
tsunami- wave dynamics. 

This analysis in the narrow tsunami case does highlight the important aspects of the dynamics 
in a clear manner, allowing a clean separation of the two cases described above. We have carried 
a full numerical integration of the equations of motion that reproduce the results described above. 
The criterion for the crossover between tsunami-wave and classical chaotic inflation is determined 
by the relative contributions to the energy density from the quantum fluctuations in the tsunami 
wave state as compared to the energy density of the zero mode . 

The previous results [eqs.(|90|)-( |10lD ] can be easily generalized for generic continuous distribu- 
tions of modes and for mixed states. One has just to make the changes indicated in eq.(|85|) for 
generic distributions and in eq.(^) for mixed states. 

The generalized slow-roll conditions takes then the form: 



f]\T)+g I q'dq\fg{T)\-' coth 



dq |/,(t)|2 coth 



during the first stage of inflation. 

The dynamics is tsunami dominated provided. 



/ 



g / q^ dq{l+q^)\fqm^ coth 



>r,2 



The effective zero mode in the second stage of inflation is now given by 

Vlffir) = r,\T)+g [ q' dq |/,(r)|2 coth ^® 



These results have been verifled by numerical integration of the full set of evolution equations 
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5 Conclusions 



In this article we have studied inflation in typical scalar field theories as a consequence of the time 
evolution of a novel quantum state. This quantum state is characterized by a vanishing expectation 
value of the scalar field, i.e, a vanishing zero mode, but a non-perturbatively large number of quanta 
in a momentum band, thus its name-tsunami-wave state. 

This state leads to a non-perturbatively large energy density which is localized in the band 
of excited quantum modes. We find that the self-consistent equations for the evolution of this 
quantum state and the scale factor lead to infiation under conditions that are the quantum analog 
of slow-roll. 

The self-consistent evolution was studied analytically and numerically in a wide range of pa- 
rameters for the shape and position of the distribution of excited quanta. The numerical results 
confirm all the features obtained from the analytic treatment. 

Under the conditions that guarantee inflation, there are two consecutive but distinct inflationary 
epochs. The first stage features a rapid fall-off of the Hubble parameter and is characterized by 
the quantum aspects of the state. During this first stage the large number of quanta in the excited 
band are redshifted and build up an effective homogeneous classical condensate. The amplitude of 
this condensate is non-perturbatively large, of 0(1/A), as a consequence of the non-perturbatively 
large number of quanta in the band of excited modes. 

The second stage is similar to the classical chaotic scenario and it is driven by the dynamics 
of this effective classical condensate, with vanishing expectation value of the scalar field. Under 
the tsunami slow-roll conditions on the quantum state, the total number of e-folds is more than 
enough to satisfy the constraints of infiationary cosmology. The band of excited wave-vectors if 
not initially outside the causal horizon, becomes superhorizon during the first inflationary stage, 
therefore these excited states do not modify the power spectrum of scalar density perturbations on 
wavelengths that are of cosmological relevance today. 

Therefore, these tsunami-wave quantum states provide a quantum fleld theoretical justification 
of chaotic (or in general large field) infiationary models and yield to a microscopic understanding 
of the emergence of classical homogeneous field configurations of large amplitude as an effective 
collective mode built from the large number of quanta in the excited band. 

In addition, we recall that it is necessary to choose an initial state that breaks the <I) ^ — <!> 
symmetry in classical chaotic scenarios ^. This is not the case here. We have infiation with 
zero expectation value of the scalar field. 

For completeness we have also studied more general states and established the important dif- 
ference between tsunami (pure or mixed) quantum states leading to inflation, and thermal mixed 
states which do not lead to inflation. 

Acknowledgements: D.B. thanks NSF for support through grants PHY-9605186, PHY-9988720 
and NSF-INT-9815064. H. J. d. V. thanks the CNRS-NSF collaboration for support. F. J. C. 
thanks the Ministerio de Educacion y Cultura (Spain) for financial support through the F. P. U. 
Programme. 

A Equations of motion in the large N limit and initial states. 

In this appendix we obtain the equations of motion in conformal time for the generalized case in 
which the initial state is determined by a density matrix. The evolution of the functional density 
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matrix is given by the Liouville equation in conformal time 



H 



9 f 

9^ 



H 



d_ 
9I' 



(102) 



where the Hamiltonian H is given by eq.(15) to leading order in the large N limit. Consistently with 
the fact that in the large N limit the Hamiltonian describes a collection of harmonic oscillators, we 
propose a Gaussian density matrix 



p[7f,*,r] =AA(T)[]exp 



TTfc • 7r_fc VTfe • 7r_fc - Bk[I ) VTfc • 7r_fc 



(103) 



The hermiticity condition = p for the density matrix impose that must be real. In addition, 
since 7f(x,T) is a real field, its Fourier components must obey the hermiticity condition 7f_fc(T) = 
7f^(T); thus, we can assume A_fc(T) = ^fc(T) without loss of generality. 

The evolution equations for Ak{T)^ ■^(^) and Bf^{T) are obtained from the Liouville eq. ( 102 ) 
where the hamiltonian is given by eq. ([l5|). We find 



Al-Bl-a\T)u;lir) 



iB^ = Bk (Ak - Al 



Mp{T) = AA,(0)e-^^ -^^S^t^^^^^-^^^^)] 



(104) 



where the prime denotes derivative with respect to conformal time T. 

The normalization factor for mixed states Mp{T) is related with the normalization factor of 
pure states Afr{T) by 

Mp{T) = Mr{T) Mriry 



where 



Mr{T) =J\fr{0) exp 



r 



dT' 



N 



Writing A^. in terms of its real and imaginary parts A^. = Ajij^+iAjj^, we find that Bf^/Aji^i;. is a 
conserved quantity. Thus, we can introduce without loss of generality the variables ^/{^^(T), Aj^k{T) 
and 6fc defined by 



Bk{T) 



AR,k{T) cothOfc 
sinh Qu 



(105) 



where 0^ is a time independent real function. 
Introducing the complex variable 



Ak = An^k + iAi^k 



(106) 
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we see that it obeys the following Ricatti equation 
This equation can be linearized defining 

Then eq. ( |107] ) implies that the mode functions fk obey 



(107) 



(108) 



u^l{T) = k' + a\T) 



6 



(109) 



where TZiT) is the Ricci scalar. 

The relation ( |108| ) defines the mode functions ^Pk(X) up to an arbitrary multiplicative constant 
that we choose such that the wronskian takes the value, 



"Pk fC -'Pkfk = '^'i' ■ 

For this choice of the Wronskian the definition (|10^ ) becomes 

1 i d 



Ak 



ln|v9fc| 



WkV 2 dT 

The mass term in eq.( |109| ) given by eq. ([T7| ) requires the self-consistent expectation value 

d^k 



(110) 



:iii) 



N 



Trp vffc • 7f_ 



1 



l_„.,(^).^l,,|.„,,(^)n2) 



2[AR,k + Bk] 2Ar,i 

Thus, the evolution equations in terms of the mode functions are given by eq. (|109|) with 



1 



A ^ f Wk 



(27r)3 a2 



coth 



Ok 
2 



:ii3) 



The evolution equation of the mode functions tpk is the same as the Heisenberg equations of motion 
for the fields, hence we can write the Heisenberg field operators as 



7r(x, T) 



d^k 



(114) 



^(27r)3 

Thus, the definition (|108|) gives the relation between Schrodinger and Heisenberg pictures, since 



the functional density matrix ( 103 ) is in the Schrodinger picture. 
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The expectation value 'ipiT') [see eq. (p^)] in conformal time obeys the following equation of 
motion P 



{T) + a\T) 



n{T) 



V(T) = 



(115) 
(116) 



Hence, the evolution equations are given by (|109|) , ( Imp and ( plSD with ( [TT^ ). 

The initial density matrix in the Schrodinguer picture is determined by specifying the initial 
values of ^_R,fc, Ai^k and 0^. We will take o(0) = 1 and parameterize the initial value of Ak as 
follows, 

^i?,fc(0)=l^fc , Ai,m=uJk{Q)6k (117) 

The corresponding initial conditions for the mode functions are obtained from eq.( |117| ) using 
eq.(108) and the Wronskian constraint eq.(llO). These are given by 



y'fc(O) 



V'fc(O) = -[wfc(O) 5k + i^k] fkiO) 



(118) 



Defining the number of particles in terms of the adiabatic eigenstates of the Hamiltonian (|15|) as 
in eq. ([27|), it is straightforward to find that the initial occupation numbers are given by 



nfc(U) - (nfc(U))p(o) p; coth 

Aujk^lk 



2 



(119) 



For any mixing parameter 0;^ 7^ the density matrix represents a mixed state since Bk ^ 0, a pure 
initial state is obtained by taking 0^, = cxd, in which case i?^ — > and the density matrix becomes 
a product of a wave functional times its complex conjugate. 

It is convenient to pass to comoving time, this is achieved by the rescaling of the fields 



iP{T{t)) = ^{t)ait) , ipkiTit)) = fk{t) ait) 
in terms of which the equations of motion are 

4>{t) + 3 H{t) 4){t) + M'^{t)>p{t) = 

fc2 



(120) 



fk{t) + 'iH{t)fk{t) + 



a^t) 



+ M\t) 



M\t) = n?+mt) + \ct>''{t) + -^ 



fk{t) = 
d?k 



(27r)2 



\fk{t)V coth 



2 



:i2i) 



where the dots denote derivative with respect to the comoving time t. The initial conditions for the 
order parameter are its initial value <^(0), and its initial derivative (/'(O). For a(0) = 1, the initial 
conditions for the fluctuations are given by 0^ and 



A(o) 



1 



A(0) = -[uJk{^)5k+H{Q)+mk] fk{0) 



(122) 



[Those are the transformed of the initial conditions in conformal time eq. (|118|). 



28 



References 

[1] M. S. Turner and J. A. Tyson, Rev. Mod. Phys. 71 (1999) S145. 

[2] See for example, M. Kamionkowski and A. Kosowsky, Ann. Rev. Nucl. Part. Sci. 49, 77 
(1999). A. H. JafFe et al., 'Cosmology from Maxima-1, Boomerang and COBE/DMR CMB 



Observations', |astro-ph / 0007333 



[3] E. W. Kolb and M. S. Turner, The Early Universe, Addison Wesley, Redwood City, C.A. 1990. 
A. H. Guth, 'Eternal Inflation', |astro-ph/0101507 . 

[4] P. Coles and F. Lucchin, Cosmology, (John Wiley, Chichester, 1995). 

[5] A. R. Liddle and D. H. Lyth, Cosmological Inflation and Large Scale Structure, (Cambridge 
University Press, 2000) 

[6] A. R. Liddle, 'The Early Universe', in 'From quantum fluctuations to cosmological structures', 
eds D Valls-Gabaud, M A Hendry, P Molaro and K Chamcham, Astronomical Society of the 
Pacific Conference Series, Vol. 126, 31 (1997). 

[7] J. Lidsey, A. Liddle, E. Kolb, E. Copeland, T. Barreiro and M. Abney, Rev. of Mod. Phys. 69 
373, (1997). 

[8] D. Boyanovsky and H J de Vega, in the Proceedings of the Vllth. Erice Chalonge School 
on Astrofundamental Physics, N. Sanchez ed., Kluwer, Series C, 2000, |astro-ph/000"6446| . D. 
Boyanovsky, H. J. de Vega and R. Holman, Phys. Rev. D49, 2769 (1994). D. Boyanovsky, D. 
Cormier, H. J. de Vega, R. Holman, A. Singh, M. Srednicki, Phys. Rev. D56, 1939 (1997). D. 
Boyanovsky, D. Cormier, H. J. de Vega and R. Holman, Phys. Rev. D55, 3373 (1997). 

[9] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman and S. P. Kumar, Phys. Rev. D57, 
2166 (1998). 

[10] M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press, 1997; D. F. 
Walls and G. J. Milburn, Quantum Optics, Springer- Verlag, NY, 1994. 

[11] L. P. Grishchuk and Y. V. Sidorov, Phys. Rev. D42, 3413 (1990) and Class. Quant. Grav. 6, 
L161, (1989). 

A. Albrecht, P. Ferreira, M. Joyce, T. Prokopec, Phys. Rev. D50, 4807 (1994). 

C. Kiefer, D. Polarski, A.A. Starobinsky, Int. J. Mod. Phys. D7, 455 (1998); 

C. Kiefer, J. Lesgourgues, D. Polarski, A. A. Starobinsky, Class. Quant. Grav. 15 L67 (1998). 

[12] R. D. Pisarski, 'Nonabelian Debye screening, tsunami waves, and worldline fermions'. Proceed- 
ings of the International School of Astrophysics 'D. Chalonge', Erice, Italy, Sept. 4-15, 1997, 
p. 195, N. Sanchez and A. Zichichi Editors, Kluwer Acad. Publ., Dordrecht, 1998. 

[13] D. Boyanovsky, H. J. de Vega, R. Holman, S. Prem Kumar and R. D. Pisarski, Phys. Rev. 
D57, 3653 (1998). 

[14] F. J. Cao, H. J. de Vega, Phys. Rev. D63, 045021 (2001). 



29 



[15] See for example, J. W. Harris and B. Muller, Ann. Rev. Nucl. Part. Sci. 46,71 (1996) and 
references therein. 



[16] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, P. R. Anderson, Phys. Rev. D50, 2848 
(1994). F. Cooper, Y. Kluger, E. Mottola, J. P. Paz, Phys. Rev. D51, 2377 (1995); F. Cooper 
and E. Mottola, Mod. Phys. Lett. A 2, 635 (1987); F. Cooper and E. Mottola, Phys. Rev. 
D36, 3114 (1987); F. Cooper, S.-Y. Pi and P. N. Standoff, Phys. Rev. D34, 3831 (1986). 



CD LO ^ CD 
+ + + CD CD CD 
O) CD CD CD CD CD CD CD 
^ -I— -I— 1— -^CO-i— to 




CD 
CD 



CD 
OO 



C3 
CD 



C3 



CD 



CD 



O3 



Figure 1: Tsunami inflation: isolines of constant number of efolds obtained from eq.(|8^) (valid for 
gT,o <C1), for m = W^Mpi, X = IQ-^^ and N = 20. 
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